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Abstract 

We set up a systematic expansion of the prepotential for M = 2 supersymmetric Yang- 
Mills theories with SU(N) gauge group in the region of strong coupling where a maximal 
number of mutually local helds become massless. In particular, we derive the hrst non¬ 
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providing a coupling between the different U(l) factors, and which governs for example 
the strength of gaugino Yukawa couplings. 


* Research supported in part by the National Science Foundation under grants PHY- 
95-31023 and DMS-95-05399, and by the Harmon Duncombe Foundation. 


1 



I. INTRODUCTION 


The leading low energy behavior of A/” = 2 snpersymmetric Yang-Mills theory is 
completely determined by the effective prepotential JF, as a fnnction of the modnli of 
ineqnivalent vacnnm states. The effective prepotential T is obtained in terms of the peri¬ 
ods of a meromorphic Abelian differential dA, on a family of anxiliary Riemann snrfaces 
parametrized by the modnli of vacna. This constrnction was pioneered in [1] for gange 
gronp SU(2), extended in [2] to SU(N), to the inclnsion of hypermnltiplets in [3] and to 
other gange gronps in [4], The power of this resnlt resides in the fact that it is valid for 
all ranges of the conpling constant, inclnding strong conpling. Of special interest is the 
dynamics of the theory close to the singnlarities in the modnli space of vacna; it is there 
that extra massless degrees of freedom appear. 

At weak conpling, the Seiberg-Witten Ansatz reprodnces the well-known pertnrbative 
expression, together with mnlti-instanton corrections. For gange gronp SU(2), an expan¬ 
sion is readily derived from [1], and was fonnd to agree with direct calcnlations in held 
theory in [5] to the orders available (np to held redehnitions). For gange gronps SU(N), 
SO(N) and Sp(N), (and any nnmber of hypermnltiplets in the fnndamental repesentation 
of the gange gronp), efficient methods for obtaining the instanton corrections to arbi¬ 
trary order were developed and 1- and 2-instanton contribntions evalnated explicitly in [6]. 
Geometrically, in the weak conpling region, the family of anxiliary Riemann snrfaces de¬ 
generates to a prodnct of two spheres joined by N thin tnbes, thns prodncing an instanton 
expansion in terms of rational fnnctions. 

At strong conpling, singnlarities in the prepotential arise along snbvarieties in the 
modnli space of vacna where certain magnetic magnetic monopoles or dyons become mass¬ 
less [1]. In view of the mass formnla for BPS states, snch massless states will appear 
whenever one (or more) of the periods of the diherential dA vanish. Geometrically, the 
strong conpling singnlarities will occnr at the snbvarieties in the modnli space of vacna 
where the anxiliary Riemann snrfaces develop nodes. 

The charge assignments of the magnetic monopoles or dyons that become massless 
at a given node, may be obtained directly from the Picard-Lefschetz theorem, and their 
associated monodromies were fonnd to pass varions physical consistency checks [1-4]. Their 
valnes may be nsed to obtain the logarithmically singnlar contribntions to the prepotential, 
associated with the renormalization gronp behavior of the infrared free light dyons at those 
points. For gange gronps of rank at least 2, Argyres and Donglas [7] fonnd that special 
snbvarieties occnr in the modnli space of vacna where mntnally non-local dyons become 
massless simnltaneonsly. 

A more systematic analysis of the prepotential aronnd strong conpling singnlarities 
has been difficnlt to obtain in general, however. One method makes nse of the Picard- 
Fnchs eqnations for the periods of the meromorphic 1-form dA. It concisely reprodnces 
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the results for gauge group SU(2), and leads to explicit expressions for the quantum order 
parameters for SU(3) in terms of Appell functions, from which a strong coupling expansion 
may be obtained [8]. For gauge groups of rank 3, the associated Picard-Fuchs equations 
were derived only recently, and their complexity and sheer size suggest that it will be 
virtually impossible to generalize the results of [8] to gauge groups of higher rank. 

In the present paper, we shall set up methods by which systematic strong coupling ex¬ 
pansions of the prepotential at strong coupling singularities may be obtained directly from 
the geometry of the Riemann surface and the meromorphic 1-form. The key ingredients are 
as follows : (1) the prepotential is easily obtained from a renormalization group equation 
derived in [9-10], (2) the ”beta-function”, i.e. the right hand side of the renormalization 
group, is obtained by residue calculations on the small cycles of the Riemann surface. 

We study theories with gauge group SU(N) and without hypermultiplets, but the gen¬ 
eralization to other gauge groups and the inclusion of hypermultiplets is analogous. We 
shall concentrate on the region in the moduli space of vacua where the maximal number 
of mutually local dyons become massless (which is N-1 for gauge group SU(N)). These are 
precisly the strong coupling singularities for which the auxiliary Riemann surface degen¬ 
erates to two spheres connected by N thin tubes, just as was the case for weak coupling, 
so that the prepotential may be expanded in a series of rational functions. These maximal 
strong coupling singularities are also the ones whose dynamics was studied in [11], in part 
in the limit of large numbers of colors. In particular, we derive the hrst non-trivial non- 
perturbative correction, which is the hrst term in the strong coupling expansion providing 
a coupling between the different U(l) factors, and which governs for example the strength 
of gaugino Yukawa couplings. 


II. EXPANSION AROUND CURVES OF MAXIMAL SINGULARITY 


The renormalized order parameters a^, their duals a^k and the prepotential are de 
hned by 

9jF(a£)) 


27riak = (p dA, 2TiiaDk = f dA 


Bi. 


dttDk 


O'k 


in terms of the meromorphic Abelian differential dA, given by 

xA'{x)dx 


dX = 


v/A(^)2^^4A^ 


( 2 . 1 ) 


( 2 . 2 ) 


Here A{x) is a polynomial of degree N in x, whose coefRcients vary with the moduli of vacua, 
and A is the renormalization scale. The associated Riemann surface (or curve) is of genus 
N—1 and is given by = A{x)‘^ . Without loss of generality (but possible rescaling 

of A), the polynomial A may be normalized so that A{x) — x^ -|- ux^~‘^ + 0{x^~^). 
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We shall make use of the renormalization group equation for the prepotential found 
in [8,9] and expressed here in terms of the dual order parameters : 


N 

(iDk 

k=i 


dJ^ 

dttDk 


- ‘IT 


2N 

27ri 


u{cId) 


(2.3) 


where u is the expansion coefhcient of A{x) as dehned above. This equation is extremely 
useful, because it now suffices to obtain the coefficient n as a function of the periods aok- 
When the Bk cycles degenerate, such a result is easily gotten from residue calculations, as 
was shown for the weak coupling regime in [10]. On the other hand, obtaining T directly 
from the the dual Ak cycles, of (2.1), would involve complicated integrals over long cycles. 
Notice that in the above equation (2.3), the roles of the renormalized order parameters 
and of their duals a^k may be interchanged. 


(a) The maximal singularity curves 

The strong coupling singularities at which a maximal number of mutually local mag¬ 
netic monopoles or dyons become massless are characterized by curves = Aq[x)‘^ — 
with a maximal number of distinct double zeros of and thus no zeros of order higher 
than 2. (Zeros of order 3 and above produce Argyres-Douglas singularities, and we shall 
not sudy those here.) As was shown in [11], these curves are given by the following (Cheby- 
chev) polynomials 

Ao{x) = 2A^Co{-^) Co{z) = cos(N arccos(z)) (2.4) 

Since only entered the original problem, and Aq(x) is a function involving only A^, 
there are in general N different solutions Aq{x), associated with the N different roots of 
unity which can multiply A. All these solutions are physically equivalent and we shall 
restrict to the one for which A is real. 

The neighborhood of the maximal singularity curve may be parametrized by the N — 1 
parameter family of polynomials P of degree N — 2, whose coefficients are assumed to be 
small in units of A. Without loss of generality, we may render its coefficients dimensionless 
by scaling out suitable factors of A, similar to the form of (2.4). We then obtain the 
following parametrization of the neighborhood of the maximal singularity 

A{x)=Ao{x) + 2A^P{^) (2.5) 

Rescaling a; by a factor of 2A so that it too becomes dimensionless, we have 


dX 


x{Cq{x) + P'{x))dx 

y(c„(i) + pw)2-i 


( 2 . 6 ) 
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Henceforth, we shall measure the renormalized order parameters and their duals 
in units of 2A and thus set A = 1/2. 

The double zeros of Co{x) are given by 

rT'TT 

Cfc = cos(—) k = !,■■■,N -1 (2.7) 

The coefficients of the polynomial P{x) parametrize the way in which the N — 1 double 
zeros are being approached. Generically, the effect of adding P{x) in (2.5) is to lift each 
double zero to a small branch cut in the double sheeted plane representation of the curve. 
We dehne the degenerating cycle H/j, k = 1, ■ • •, A^ — 1, as the cycle that surrounds the 
branch cut created from the double zero Cfc once in the counterclockwise direction. The 
cycles Ak can then be chosen to be any set of conjugate cyles, but we shall not need them 
here. 


(b) Expansion aronnd degenerating cycles 

We assume that the coefficients of P{x) are small compared to 1, and that the curves 
Bk are at a distance from Cfc of order a fraction of 1, (and thus much larger than the 
coefficients of P{x)). Then, we may expand the period integrals for the a^k periods 
associated with the degenerating cycles in a convergent power series of P{x). The 
radius of convergence for this expansion can be shown to be of order a fraction of 1 in 
these units. 

To perform this expansion in an efficient way, it is most convenient to introduce an 
auxiliary parameter fi which interpolates between Cq and Cq-\-P by C^{x) = Co{x)+iJiP{x). 
The period integrals now become 


0,Dk = — 



dx P{x) 

^ ^/C^ixy -1 


oo 



m=0 


( 2 . 8 ) 


with the expansion terms a 


(m) 

Dk 


given by 


a 


(m) 

Dk 


^ ’ r(i)m! 



dx P{x)‘^~^^{lj.P{x) +2Co{x))‘^ 
27rf (C'o(a:)2 - l)^+h 


(2.9) 


We £x the square root sign by {Co{xY — f sin(A^arccosa:). 

For low rank A^ — 1 of the gauge group, the polynomial P{x) will be of low degree, and 
one may use the coefficients of the polynomial as suitable classical parameters of moduli 
space of vacua. As soon as N becomes larger, this parametrization is not convenient 
however. To identify a good parametrization, we analyze the hrst term (m = 0) in the 
above expansion. We introduce the very convenient notation 


Pk = P{Ck)\ 


^k 


sin( 


kn 


Ck = cos( 


kn 


k = l,---,N-I (2.10) 
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The first expansion term is given by 


“dI = Pt ( 2 - 11 ) 

which indicates that convenient local coordinates of the modnli space of vacna aronnd the 
maximal singnlarity are the valnes Pk of P{x) at the donble points. We recover P{x) with 
the help of the Lagrange interpolating formnla giving the polynomial of degree — 2 in 
terms of its — 1 valnes Pfc at A^ — 1 hxed points Ck 




fc=i 


(x - Ck)Cn[ck) 


( 2 , 12 ) 


To obtain nsefnl formnlas for the expansion coefficients, we take advantage of the fact that 
the factor of Cq{x) almost eliminates the denominator in the integrand of (2.9). More 
precisely, we have 


Cq{x) iN 

Vl - x‘^ 


(2.13) 


We now introduce the functions 11^(2:) and Gk{x), which are analytic at x = Cfc; we also 
dehne the functions n(a:) and G{x), which have simple poles at all Ck, and which we shall 
make use of later on. 


Uk{x) == (x - Cfc)n( 2 :) 


Gk{x) = (x - Ck)G{x) 


n(a:) = 


P{x) 


G{x) = 


Goix) 

G(,{x) 


N-l 


-EH' 


X - Cl 


N-l 

xN — 

1=1 



(2.14) 


In terms of these functions, the expansion coefficients of (2.9) are easily evaluated, 
and we obtain 


a 


(m) _ 
Dk ~ 


. r(m + i)/r(i) 

t - - - - — 

771/^ ^ 



dx (1 — x'^) 2 

27rf {x — 


Ukixr+^ 


ifiUk{x) + 2Gkix)r 


Um+iyni) f' 

m!(2m)!iV2’"+i 


(1 — 

X=Ct, 


--^Uk{xr+\pUk{x) + 2Gk{x)r 


(2.15) 

Here, the x-integral was evaluated by residue methods. Formula (2.15) provides an ex¬ 
pansion to all orders in Pk of the dual quantum order parameters aok- To evaluate the 
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expansion coefficients, it will prove convenient to have the following derivative terms 


Giick) 

G'fc(cfc) ={-N'^sI + -si + 2 ^ 1 )^ 
G'kick) ={-N‘^sI + -cl + 
as well as the fact that Ilk{ck) = —{—)^slPk- 


(2.16) 


(c) Contributions to second order 

We shall now use the above expression to evaluate contributions to second order in 
Pk hrst, and hnd that 


Dk quadratic in Pj. 


■^py-G-4-lP^N^4)PS-Kic.r- 


+ i-fslPtnUct) + 3(-)''ctPtni(ct)] 

The derivatives are easily evaluated and given as follows 

N-l 2 

pk « 


N-l 

K.(ck)=2j2(-)‘ 

l^k 


(Ck - Cl)' 


:Pl 


(2.17) 


(2.18) 


Higher order terms may be calculated analogously, but we shall not need them here. 
Instead, we shall use some shortcuts and compute directly the beta function. 


III. THE BETA FUNCTION AND THE PREPOTENTIAL 

The renormalization group beta function u, discussed in (2.3) is obtained in terms of 
the leading coefficient multiplying x^~‘^ of P{x). Remarkably, this quantity also assumes 
a simple expression in terms of the coordinates Pk, and we hnd 

N-l 

«=-4A=‘fi; u^-J^i-fslPk (3.1) 

k=l 

Here, we have dropped a constant contribution (—AA^) to u, since it would yield just a 
constant shift of the prepotential P which is physically immaterial. 
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Now, instead of calculating out the expansions of in terms of the classical variables 
Pk, then inverting to obtain an expansion of Pk in terms of a^fc, and hnally inserting the 
expansions of Pk into (2.19) to obtain the beta function u, it is much more advantageous 
to exhibit the summation that occurs in u before effecting the substitution. In particular, 
this procedure will avoid having to perform major resummations which would have to be 
carried out by contour integrations anyway. Thus, we shall instead consider the following 
sum which directly exhibits u 


N-l 


U 


f ^ ^ ^kOtDk T ^ ^ 


a 


(m) 


D 


a 


(m) 

D 


N-l 


E (m) 

^ka^Dk 


(3.2) 


fc=l 


m=l 


k=l 


The resummed expansion parameters may be expressed in terms of the functions 

n(a:) and G{x)^ dehned in (2.14), and are given by 


(m) 


N-l 


an — 


r(m + i) 


^ r(i)m!iV2m+i 


dni/ 


dx 


^k 


- Il{x)^+^{^ill{x) +2G{x)Y 


Bu 


21X1 (l_a;2)m+2 


(3.3) 

What prevents us from carrying out the sum over k by connecting the various contours 
Bk in the above formula is twofold. First, the square root would force us to perform a line 
integration around the lines from ±1 to cx). This would not represent an advantage over 
a direct evaluation with residue methods. Second, the factor Sk depends upon /c, so we 
cannot just add up the various contours Bk- 


(a) Combining contours 

We can solve both problems by replacing Sk inside the above integrand by a k- 
independent function that takes the value Sk at x = Ck for all /c = 1,---,A^ — 1, but 
whose derivatives of orders 1, ■ ■ ■, 2m vanish at this point. Replacing Sk by such a function 
inside (3.3) for each value of m will not alter the value of the integral. It is very easy to 
obtain such a f un ction by expanding the value at Ck in a Taylor series around x instead. 
For example, dehning s((r) = Vl — we have 


Sk = s{Ck) = s{x) + (Cfc - x)s'{x) + ^(Cfc - xYs"{x) + 0{{Ck - xY) 


To order m = 1, the expressions we need are given by 
Sk 1 


(1 — x’^Y'G 1 — x" 


+ (^ Cfc) (1_ ^2)2 2^"^ {\-x‘^f 


+ 0((i - Ct)®) 


(3.4) 

(3.5) 


while to order m = 2, we have 


(1 — 


G{x) 


X 


(1 


2\2 


{1 — X 
2 


2'|2 


G{x) - - 


1 1 


2 1 


(3 + llCfc + 28iV^Cfc)(a: - CkY + 0((x - CkY) 


x^ 

2 


^(l + 2N^)(x-Ck) 


244 


(3.6) 



Inside the integrals, the 0{) terms vanish identically and may be dropped. What we have 
gained is that the hrst term in (3.5) and the hrst three terms in (3.6) are /c-independent 
and varions contonrs in the snm (3.3) may now be collected together. As can be established 
by inspecting (3.3), (3.5) and (3.6), the only poles other than those occnring at c^, k = 
1, ■ ■ ■, N — 1 are the poles at a: = ±1. Thns, we have for those parts of the integral 



(3.7) 


The remaining terms are handled by standard residne calcnlations. 


(b) Contributions for m = 1 

Substituting (3.5) into (3.3) for m = 1, we hnd, using also (3.7) 


.,(1) 


1 


dx 1 


1 . 


^ M +n(x)^G(x)) 


+ 


iV-l 


2N^ ^ Jn 2TTi 
k=i 


dx 2n 1/ M + ^k{xYGk{x) 

-Ml- 


(3.8) 


It is straightforward to carry out the residue calculations at ±1 in the hrst integrals, and 
at Cfc, /c = 1, • • •, A^ — 1 in the second integrals. Using also the fact that G(±l) = ±1, we 
hnd 


1 r 1 I 

- IVs E - 4(-)'''ciP,nUci) + (-)'■■(- + 2cl)Pl - 2cMXM) 

k=l ^ ^ 


Notice that the combinations n(±l) are simply given by 

N-1 

(3.9) 

k=l 

(3.10) 


Making use of expression (2.18) for n'^(cfc), and some straightforward trigonometric 
identities, we hnd 


a 


( 1 ) 

D 


N-1 


4iV3 


^E4A +Y^(n(i)^-n(-i) 


k=l 
N-1 


' E + 24 )P| - 2cMKM} 


4iV3 


k=i 


(3.11) 
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Notice that the quadratic part is diagonal in P^. 


(c) Contributions for m = 2 

In a manner analogous to the computation of the m = 1 contributions, we get those for 
m = 2 by substituting (3.6) into (3.3) for m = 2. Using the fact that G'(±l) = |(2A^^ + 1), 
we hnd 


( 2 ) 

ay = 


1 


24iV5 
N-l 

E 

k=i 


{2N^+ 


1 




+ 


12 


-N^cDPi 


(3.12) 


(d) Expressing P^ and the beta function in terms of aDk 

To obtain the beta function to cubic order in it suffices to substitute into (3.11) 
and (3.12) the expansion of Pk in terms of aok up to quadratic order only. This inversion 
is easily obtained from (2.11) and (2.17), and we hnd 


Pk — —{ — ) O'Dk 

^k 




0 2 


2 + ];^nfc(cfc)^ 


^3 

T *^UDfcn^(cfc) + —SfcCfca£)fcn^(cfc)] + 0{a£)k) 


(3.13) 


Putting all together, we hnd for the beta function the following expression 


N-l N-l 

a = - i E - iv E - ^(n(i)= - n(-i)“) 

k=l k=l 


(3,14) 

Here, the derivatives of H^ are evaluated as follows 


N-l 


Ei 




/ 7 2^ 

<I2 




2i 

'^'ok^ki'^k) — -rrCkSkoj^i^Iiyck) 


N 


nyck) ^iN 

K{ck) = - 


SlClDl 


2iNj2 


l^k 


siaDi 

{ck - Cl) 


2 


(3.15) 


It remains to rearrange the sum in the last line of (3.14), and to combine with the n(+l)^ — 
n(—1)^ term of the hrst line. The result is as follows 


N-l 

U i ^ ^ ^kC^Dk 
k=l 


1 

4iV 


N-l 


Y1 ^Dk 

k=l 


i 

16iV2 


pAT-l 

E 


*■ k=l 



N-l 


^E 


(^‘pk^pPi 
(ck - q)^. 


(3.16) 
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Restoring the original A-dependence, we now easily find the renormalization gronp beta 
fnnction 


N-l 


N-1 


rN-1 


U = 


2 *A '^2 SkClDk + '^2 ^‘ok + 


k=l 


k=l 


32N‘^A 


N-l 




k=i 


k^l 


(Ck -Ci)^ 


(3,17) 


(e) The prepotential 

Starting from the resnlt for the beta f un ction in (3.17), and making nse of (2.3), we 
obtain onr final expression for the prepotential at strong conpling 


jF(aD) 


2NA 

7T 


AT-l 

skttDk 

k=i 


i 

47r 


+ 


1 

327riVA 


pW-l 

E 


^ fc=i 


N-l 

E 2 1 O'Dk 

k=l '' 



N-l 


^E 


OpkO'DlSl 
{Ck - q)2_ 


(3.18) 


Here, the valnes of the normalizations Afc are not available directly from the integration of 
the renormalization gronp eqnation in (2.3). (In principle, the renormalization gronp eqna- 
tion will determine the prepotential T throngh (2.3) only np to a homogeneons fnnction 
of degree 2 in Here, however, we are performing an expansion in powers of a^jt/A, 
and we know what that expansion looks like from the evalnation of the A/j periods. In 
particnlar, only a single logarithmic term will appear. Thns, knowing the form of the ex¬ 
pansions of both u and JF, the only remaining indeterminacy of the renormalization gronp 
eqnation are terms qnadratic in aDk-) They can be obtained from a direct calcnlation to 
that order of the Ak periods, a calcnlation that has already been carried ont by Donglas 
and Shenker [11]. We shall qnote their resnlts here : 

log^ = ^+logSfc (3.19) 


(f) Comparison with known results 

For SU(2), N=2, we have si = 1 and thns 


_ 4A i 2 , I-? 


(3.20) 


which agrees with [1,8] after snitably adjnsting for the slightly different definitions of A. 
For SU(3), N=3, we have si = S 2 = v^/2 and ci = —C 2 = 1/2 and thns 


JF =- ^ (am -I- aD2) 


* / 2 1 ^D1 . 2 1 _ ^D2 


TT 


A, 


+ a^2log 


A. 


+ 


1 


(3.21) 


144V37rA 


{a-Di + «D2)(4a|)]^ — l?>aDiaD2 + 4a^2) 
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Again, taking into account the difference in definition of A with [8], we hnd agreement. 
Finally, for all SU(N), the logarithmic terms were computed in [11], and the normalization 
of the logs agrees. 

(g) Physical interpretation 

To understand more clearly the physical origin of the non-perturbative correction 
terms that we have calculated, it is helpful to derive, within this approximation, the (dual 
gauge) coupling constant /daomdaDn- We hnd 


D 


r. 


r" rn.n. 


^ 1 O^Dm 

—log —— 
2 tI Km. 


O^Dm 


3z 3 
47r IGttA^A 


E 


2aDiSi 


SttNA (Cm 


ciY 


1 T (^Dn^r 


(3.22) 


47riVA 


[ Crrj, Cr 


Identifying r^m terms of the dual coupling constant g^m and dual instanton angle 9Dm, 
as usual, we hnd that the cubic correction terms in the prepotential have a dependence on 
gom and Oom through 


^27ViT 


D 

mm 


— e 


- \-iODm 


(3.23) 


suggesting that the ehects are due to (magnetic sector) instantons. It would be very illu¬ 
minating to see how the same ehects can be identihed from the physics of light monopoles 
near the strong coupling points. 
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